Complex dynamics in simple systems with seasonal parameter oscillations 
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We study systems with periodically oscillating parameters that can give way to complex periodic or 
non periodic orbits. Performing the long time limit, we can define ergodic averages such as Lyapunov 
exponents, where a negative maximal Lyapunov exponent corresponds to a stable periodic orbit. By 
this, extremely complicated periodic orbits composed of contracting and expanding phases appear 
in a natural way. Employing the technique of e-uncertain points, we find that values of the control 
parameters supporting such periodic motion are densely embedded in a set of values for which the 
motion is chaotic. When a tiny amount of noise is coupled to the system, dynamics with positive 
and with negative non-trivial Lyapunov exponents are indistinguishable. We discuss two physical 
systems, an oscillatory fiow inside a duct and a dripping faucet with variable water supply, where 
such a mechanism seems to be responsible for a complicated alternation of laminar and turbulent 
phases. 

PACS numbers: 05.45.Ac,47.60.+i 
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I. INTRODUCTION 

When studying deterministic dynamical systems, it 
has become practice to distinguish between chaotic and 
(quasi-)periodic solutions, where chaos has been seen as 
a novel and strange behavior. We study a class of dy- 
namical systems, where such a distinction is not useful. 
Depending on control parameters, our systems will either 
generate stable periodic orbits or they will behave chaoti- 
cally, but these two types of motion will be essentially in- 
distinguishable from each other both on computers (due 
to finite precision) and in real experiments (due to weak 
noises). 

Seasonal variations are a prominent source of slow peri- 
odic parameter fluctuations in biological, ecological, geo- 
chemical and geophysical systems. But also in many 
technical and physical situations, slow periodic oscilla- 
tions of system parameters do occur. Speaking of time 
dependent parameters, we imply that there is no feed- 
back from the system under study to the variation of the 
parameters, whose time dependence can either be con- 
sidered as given (non-autonomous situation) or can be 
ruled by its own periodic autonomous dynamics. More- 
over, we focus on situations where the typical time scales 
related to the system dynamics for fixed parameters is 
much faster than the time scale related to the parameter 
variation, as it is typical of many processes subject to 
seasonal variations. The opposite case, where both time 
scales are comparable _has been studied as an open loop 
control mechanism 111 



In what follows we discuss a scenario, where slow har- 
monic parameter variations introduce an alternation of 
expanding and contracting phases. We will show that er- 
godic averages can be performed as usual and hence the 



motion in the long time limit is clearly classified to be 
either periodic or chaotic. However, as our analysis will 
show, in practical applications, chaotic motion will be 
indistinguishable from periodic motion, both in numer- 
ical simulations and in experiments. In the latter case, 
this is another example of "stable chaos" Q. This indis- 
tinguishablility implies a robustness of the phenomenon 
despite the fact that there are parameter regimes where 
stable periodic motion and chaos are both supported by 
a dense set of parameters. 

This rather unexpected behavior is shown to exist in 
numerical experiments of an oscillating flow inside a duct 
and a dripping faucet with variable liquid supply. 



II. MAP WITH PERIODIC PARAMETER 
VARIATIONS 

Many physical experiments and the corresponding 
model systems possess solutions which, due to dissipa- 
tion, relax to rather uninteresting fixed points. Such sys- 
tems are often exposed to a periodic driving, such as 
electric resonance circuits. In particular when the sys- 
tem without driving has a two-dimensional phase space, 
chaos can only appear with the driving term. 

This paper deals with a very different class of driven 
systems; we assume that our system without driving can 
behave chaotically, depending on the values of a control 
parameter. This parameter is then varied periodically, 
with a period which is much longer than the time scale of 
the internal dynamics on which the auto-correlation func- 
tion in the chaotic regime would relax, or which would 
govern the relaxation to a stable fixed point in parameter 
regimes where it exists. Hence, the temporal dependence 
of the control parameter has no influence on the short 
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time dynamics of the system, but it causes transitions 
between different dynamical behaviors which might exist 
for different values of the control parameter. 

In what follows we choose the quartic map / given by 

x,i+i = f{a,Xn) = 4aa;„(l - x„)[l - axn{l - Xn)] (I) 

with < a < 4. This map has the interesting property 
that for a = flc = 3.375 a tangent bifurcation occurs 
where out of a chaotic region that covers the whole [0, 1] 
interval a period one window appears. However, other 
maps can be also be considered, for example the logistic 
map for values of the control parameter near a period 3 
window. 

Imposing a temporal dependence on the control pa- 
rameter leads us to consider the two dimensional skew 
autonomous system 

Xn+i = /(a(0„),x„) (2) 
= 4a(<^„)a;„(l - a;„)[l - a(0„)a;„(l - a;„)] 

where a{(j)n) = ao+ai cos 0„, ag is chosen near Uc, ai ^ 1 
and Lu — 2tt/N, N ^ 1. A small value of uj leads to a 
well defined time scale separation between the parameter 
variation and the dynamics of the system while the choice 
of flQ and fli ensures that a oscillates from values where 
the unperturbed map of Eq. ^ is chaotic to values where 
it has a stable fixed point. 

A typical sequence of iterates of Eq. Q (after discard- 
ing a transient) is shown in the top panel of Fig. ^ We 
see an alternation between irregular fluctuations of x and 
regular episodes due to the oscillating values of a shown 
in the bottom panel. Due to the slowness of the change of 
a, the trajectory can relax toward the stable fixed point 
for those values of a„ where it exists, whereas it follows 
an irregular trajectory for a„ < 3.375, where the map of 
Eq. ^ is chaotic. 

Since the (/)-dynamics is not mixing, this system 
clearly has (at least) one invariant measure in its two- 
dimensional phase space for every (pQ, on which ergodic 
averages are well defined. The Jacobi matrix of Eq.Q 
is triangular, which immediately shows that one Lya- 
punov exponent (corresponding to the 0-dynamics) is 
zero, whereas the other A = (A„) is found as an aver- 
age over an infinitely long trajectory of the instantaneous 
stretching factors A„ (middle panel of Fig. ^ defined by 

\ a/(a((/)„),a;„) 
\i — iog . 

OXn 

Assuming ergodicity, A can also be found as an aver- 
age of the stretching factor over the invariant measure. 
The seemingly intermittent dynamics shown in in the 
upper panel of Fig. ^ has a well defined non-trivial Lya- 
punov exponent, whose sign depends on the details of 
the a;-dynamics. The alternation between unstable and 
stable phases is different from the typical intermittency 
scenarios The contributions of the stable and un- 
stable phases to the Lyapunov exponent have a sensitive 
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FIG. 1; Two hundred iterates of the map Eq. ||5J for ao = 
3.4, ai — 0.1, uj — 27r/100 (top panel), the instantaneous 
stretching factors A„ (middle panel) and the parameter a„ 
(bottom panel). 



and subtle dependence on the choice of oq and ai as we 
show in Sec. IVII 



III. INDISTINGUISHABILITY BETWEEN 
CHAOTIC AND NON CHAOTIC ORBITS 

For the map of Eq. ^ almost any initial condition with 
a > ac, \a — ad 1 settles after a transient on a fixed 
point i(a). On the other hand, for a < ac, \a — ad 1, 
the invariant set is one-dimensional for a dense set of pa- 
rameter values, but shows type-I intermittency because 
of the closeness to the tangent bifurcations. For the map 
of Eq. ^ during one period of the auxiliary variable </>„ , 
the parameter a{(j)n) alternates between the periodic and 
chaotic regime of the map of Eq. (QJ and therefore A 
has contributions with negative and positive signs. De- 
pending on which contribution has a larger modulus, the 
overall dynamics is either chaotic or not. This in turn, 
depends on the values of ao, ai and lo. 

During the iterates where the trajectory looks regu- 
lar and is near to a; (a), the fixed point of the map of 
Eq. (^1, the tangent space dynamics is contracting, and 
we call the accumulated contraction factor /c. During 
the iterations when the trajectory looks irregular, its tan- 
gent space dynamics is essentially expanding, and we call 
the accumulated expansion factor /g. Then, roughly, the 
Lyapunov exponent is A ~ (In j/d -I- In |/e|), where (■ • • ) 
now denotes the average over successive periods of (j). If 
we start two trajectories with a distance e at the begin- 
ning of the irregular phase, at its end their distance is 
e/e- This distance will shrink during the regular phase. 
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and at its end will be efefc- Hence, if /e/c on average 
is smaller than unity, two trajectories will approach each 
other, and finally will be indistinguishable. 

This alternation between contracting and expanding 
episodes is clearly visible in Fig. ^ A stable periodic 
orbit thus is as irregular as a chaotic solution, but its 
irregular segment repeats itself exactly in every period of 
oscillation of the parameter a, whereas for a trajectory 
with a positive Lyapunov exponent it does not. Hence, 
this system can create arbitrarily complicated periodic 
orbits, since by the choice of uj one can determine the 
period length and also, how many points of the periodic 
orbit are in the irregular regime. When the (/>„ dynam- 
ics is quasi-periodic instead of periodic (by choosing A'^ 
as an irrational number), also orbits with negative Lya- 
punov exponent have non-repeating irregular segments. 
By visual inspection, these cannot be distinguished from 
orbits with a positive Lyapunov exponent. 

IV. THE STROBOSCOPIC VIEW 

The special dynamics of (/)„ hinders the 2-dimensional 
map of Eq. (|2J) of having a fixed point. The shortest peri- 
odic orbit can have length N when lo = 2'k/N . Therefore, 
it makes sense to study the composition of N successive 
iterates F{x) = Q)n=if{0'{<l>n),Xn)- When the Lyapunov 
exponent of Eq. Q is negative, F should have a sta- 
ble periodic orbit or fixed point, whereas it has only un- 
stable periodic orbits and chaotic solutions for positive 
exponents. In Fig. |21 we show the graph of F for sev- 
eral values of oq. For a > Oc, a — Oc <C 1 F{x) has an 
almost super-stable fixed point. These fixed points are 
generically born and eliminated by tangent bifurcations, 
together with their unstable counterparts. For every ini- 
tial phase 00 we have a different stroboscopic map F{x). 
However they are all topologically conjugate. For our 
choice of the variation of a{4)n), however, the existence 
or absence of a stable fixed point can best be seen when 
4>o « 0, where each orbit of the system Eq. ^ assumes 
values close to the fixed point of f{x,a) in Eq. ^ for 
fixed a « flQ -t- fli . 



V. NOISE AND ROUND-OFF EFFECTS 

The results described above, and in particular the dis- 
tinction between motion corresponding to negative and 
positive Lyapunov exponents, is correct only in the ab- 
stract mathematical setting. On a computer, the finite 
precision of the internal representation of real numbers 
can enforce the motion onto a complex periodic orbit 
although its Lyapunov exponent is positive. This hap- 
pens, when the contraction factor fc during a contract- 
ing phase is too strong, so that at the end of this phase 
the trajectory has no memory of the previous expanding 
phase. For typical values of and ai we found nu- 
merically that orbits with a positive Lyapunov exponent 
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FIG. 2: The graph F as defined in the text for ao — ac = 3.375 
(full curve) showing a tangent bifurction, for ao — ac — 0.0005 
(dashed curve) showing many unstable fixed points and for 
ao — ac + 0.0005 (dotted curve) showing an almost super 
stable fixed point. 



became periodic when the contracting phase contained 
more than thirty iterates. If these phases are shorter, ei- 
ther because the period A^ is small enough, or because og 
and ai are chosen to be inside the chaotic regime, orbits 
with positive Lyapunov exponents are non-periodic as 
expected. Hence, with finite precision and large A^, one 
cannot decide, without computation of the Lyapunov ex- 
ponent, whether the system has a stable periodic orbit 
or not. 

In an experimental realization, instead of computer 
round-off errors there is external noise coupled into the 
system. Eq. (QJ has then to be modified by adding white 
noise cr^„, where < cr < 1, (Cn^n') = Sn,n' and (^„) = 0. 
This has no visible effect on chaotic solutions of Eq. , 
but it does destroy the periodicity of stable periodic solu- 
tions. Inside the expanding and hence irregular sections, 
noise is exponentially amplified and creates orbits which 
appear chaotic. Systems like Eq. ||2Jl therefore have pe- 
riodic orbits which are extremely sensitive to external 
noise, despite linear stability. 



VI. PARAMETER DEPENDENCE 

The rather complex sequence of tangent bifurcations 
leading to the creation and destruction of the stable pe- 
riodic orbits causes the orbits to depend sensitively on 
the system parameters, so that there is a complicated 
hipping from periodic to chaotic motion as a function of 
every single parameter as illustrated in Figs.OHHl In the 
first one, we show the Lyapunov exponent A as a function 
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FIG. 3: The Lyapunov exponent A as a function of the pa- 
rameter ao with ai = 0.1, A*' = 100. The inset shows a smaller 
interval of ao at a higher resolution. 



of ao (fli and N fixed). The rapid change from chaotic 
to stable periodic solutions leads us to assume that both 
types of behaviors are supported by a dense set of pa- 
rameters. The same follows from Fig. 0] where we show 
the Lyapunov diagram where white (black) dots corre- 
spond to a positive (negative) value of A There are 
clearly two regions where A has a definite sign and an 
intermediate one where slight changes of ag or a\ have 
a dramatic effect on the dynamics. We determined the 
fractal dimension of the boundary between stable and 
unstable solutions by finding the scaling of the number 
of e-uncertain points as e is varied A point (oq, oi) 

is e-uncertain if in a neighorhood of radius e there exists 
at least one point where the Lyapunov exponent has an 
opposite sign to the one evaluated at (ao,ai). In Fig.|Sl 
we show those points of the previous figure that are e- 
uncertain for e = 10~^. From the scaling of the number 
of uncertain points with e we found the box counting di- 
mension Do of the boundary between chaotic and stable 
solutions to be ~ 1.985. This result gives a quanti- 
tative support to the idea that both stable and unstable 
orbits have a dense set of parameters. We also found that 
Do ~ 2 for iV > 100 and also for iV < 20 while it has a 
minumum value Dq ~ 1.8 for TV 30. 



VII. OSCILLATORY FLOW IN A DUCT 

As a first example that displays the behavior discussed 
previously we briefly present results of numerical exper- 
iments of an oscillatory flow in a duct filled with fluid. 
This flow can be generated by imposing oscillatory pres- 
sure or velocity fields at the ends of the duct with suitably 
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FIG. 4: Lyapunov diagram with A'^ = 100. Black dots corre- 
spond to a negative Lyapunov exponent, white to a positive 
one. A total of 65,536 dots are drawn. A cut at a\ =0.1 
corresponds to the previous Fig. 
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FIG. 5: Black dots represent the e-uncertain points with 
e = 10"^ and A = 100. 



defined phase lags. The stability of these fiows can be de- 
scribed in terms of two nondimensional parameters, the 
oscillatory Reynolds number i?5 and the Stokes param- 
eter A. These are defined by Rs = U5/v and X — D/5 
where U is a characteristic velocity, v the kinematic vis- 
cosity of the fiuid, D a characterisic distance of the duct 
and S the Stokes penetration depth. This last quantity 
is defined hy 6 = ^ v /2lo with uj the frequency of the 
oscillation. It is a well established fact that zones of dis- 
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FIG. 6; Velocities u and i; as a function of time t in arbitrary 
units at the center of a two dimensional duct with expan- 
sions at the ends and with an aspect ratio (total length/cross 
section) of 20, Rs = 1521 and A ~ 4. 



tinct dynamical qualitative behavior can be identified in 
the {Rs,X) space. Specifically, it has been experimen- 
tally observed that for A > 2 and Rs < 500, the flow is 
laminar, while for Rs > 500, the flow inside the duct is 
laminar for the phase intervals where the velocity is small 
while bursts with a frequency much larger than the forc- 
ing appear near the end of the acceleration phase 8] . As 
Rs is increased, the phase interval where high frequency 
oscillations are present gets larger, but it never covers 
the cycle entirely. The origin of the high frequency oscil- 
lations is the generation of vortices due to the instability 
of the laminar flow. Numerical results agree with the 
experiment j^. 

In particular. Fig. shows the axial u and transversal 
V velocities in the middle of a duct. As can be observed, 
when the velocity is close to zero, the trace is smooth, 
indicating laminar flow, however, high frequency oscilla- 
tions appear when the velocity reaches its maximum ab- 
solute value in each cycle. Taking ujy = U/D, as the lower 
limit of the characteristic frequency of the vortices, the 
ratio of the vortices frequency to the forcing frequency, 
LUy/uj is 2Rs/X which for this example is 768. This indi- 
cates that the changes in the driving force are slow com- 
pared with the internal vortex dynamics corresponding 
with the conditions discussed in the previous sections. 



VIII. DRIPPING FAUCET WITH VARIABLE 
SUPPLY 



Our second example describes the dynamics of a drip- 
ping faucet with variable liquid supply. This problem, 



with constant liquid supply has been studied exten- 
sively ^3 • The time interval T between successive drops 
shows a complicated bifurcation diagram as the water 
supply e increases. There are two studies that are impor- 
tant in the context of the present analysis. The model 
presented by Fuchikami et al [Tll | is relevant because it 
is built on sound physical phenomenology, but unfortu- 
nately is not simple and long time calculations involving 
many drops are extremely computing intensive. On the 
other hand, the model presented by CouUet et al 
which is based on the former, is useful since it is simple 
and can be used for exploring the long time behavior. 
Both models can be adapted to analyze the system when 
the liquid supply varies as a harmonic function of time 
and they display a similar behavior with constant and 
variable liquid supply. 

The CouUet model with constant liquid supply e has a 
complicated bifurcation diagram and we choose a value 
Ec in such a way that to its left there is are period one 
solutions, and to its right a densely chaotic region. We 
now consider a variable liquid supply that varies har- 
monically around a value near Cc with a period which is 
approximately 1000 times larger than the characteristic 
period of the subsequent drop release with constant liquid 
supply and with an amplitude Ae chosen in a way that 
the water supply does not extend outside the period four 
window. These conditions allow the system to visit alter- 
natively a zone of irregular behavior and a zone of period 
four in the map of constant liquid supply. The result of 
the periodic variation of the water supply are shown in 
Fig. 13 The continuous sinusoidal line represents the liq- 
uid supply and the dots the values of the time intervals 
between successive drops. As can be clearly seen, the 
system presents zones of stable and unsatble behavior as 
for the map of Eq. (O . 



IX. CONCLUSIONS 

We discussed systems with periodic parameter fluctua- 
tions which are driven from regular to chaotic motion and 
back. Although simple in its construction, this type of 
dynamics creates very complicated orbits with a complex 
dependence on control parameters. In computer simula- 
tions and in real experiments, it is impossible to distin- 
guish the existence of stable complex periodic solutions 
and of chaotic solutions in the underlying model, since 
round-off errors and noise interact with the dynamics. 

We presented two numerical experiments that illus- 
trate the change of behavior due to a periodic variation 
of a parameter with a time scale much larger than the 
natural time scale of the system where the discussion of 
the impossibility of distinguishing chaos from order is 
relevant. 

Our detailed analysis was based on maps, but all fea- 
tures are found as well in systems with continuous time. 
It is also not relevant to assume sinusoidal variation of 
the parameter, so that we expect such a behavior to be 
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rather widespread. 
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